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ABSTRACT 

We use the shear catalog from the CFHT Stripe-82 Survey to measure the subhalo 
masses of satellite galaxies in redMaPPer clusters. Assuming a Chabrier Initial Mass 
Function (IMF) and a truncated NFW model for the subhalo mass distribution, we 
find that the sub-halo mass to galaxy stellar mass ratio increases as a function of pro¬ 
jected halo-centric radius from Msub/Mgtar = 4.43^2^23 ‘^p ^ [0.1, 0.3] h“^Mpc 

to Msub/Mstar = 75.40^19.09 '^p ^ [0-6, 0.9] h“^Mpc. We also investigate the depen¬ 

dence of subhalo masses on stellar mass by splitting satellite galaxies into two stellar 
mass bins: 10 < log(Mstar/h“^M 0 ) < 10.5 and 11 < log(Mstar/h“^M 0 ) < 12. The 
best-fit subhalo mass of the more massive satellite galaxy bin is larger than that of 
the less massive satellites: log(Msub/h~^M 0 ) = 11.14^0^73 (-^sub/Mgtar = 
versus \og{M,^b/h-^MQ) = 12 M+°aI (Vub/M,tar = 


1 INTRODUCTION 


In a cold dark matter (CDM) universe, dark matter haloes 
form hierarchically through accretion and merging(for a re¬ 
cent review, see |Frenk fc White] ( |201^ ). Many rigorous and 
reliable predictions for the halo mass function and subhalo 
mass functions in CDM are provided by numerical simula- 


tions(e.g. |Springel et al.||2008| |Gao et al.| 

201I| |Colm et al.| 

20001 Hellwing et al.|2015 Bose et al. 2016 

1. When the small 


haloes merge to the larger systems, they become the sub¬ 
haloes and suffer from environmental effects such as tidal 
stripping, tidal heating and dynamical friction that tend to 


remove the mass from them and even disrupt them (iTor- 

men et al.||1998| |Taffoni et al.||2003| |Diemand et al. 

2007 

Hayashi et al. 20031 Gao et al.||20041 Springel et al. 

2008 

Xie & Gao 

2015 

1 . At the mean time, the satellite galax- 


ies that reside in the subhaloes also experience environmen¬ 
tal effects. The tidal stripping and ram-pressure can remove 


the hot gas halo from satellite galaxies which in turn cuts off 
their supply of cold gas and quenches star formation) Balogh| 


et al.|2000l|Kawata fc Mulchaey|2008[[McCarthy et al.|2008| 

Wang et al.|2007 Guo et al.|2011 Webel et al. 12013 1. Satel¬ 


lite galaxies in some cases also experience a mass loss in the 
cold gas component and stellar component during the in- 


et al. 1999 

IGhung et al. 20091 

Mayer et al.||2001| Klimen- 

towski et a 

.|2007 IKang & van c 

en Bosch|2008 Chang et al. 


20131. 


Overall, the subhaloes preferentially lose their dark 
mass rather than the luminous mass, because the mass dis¬ 
tribution of satellite galaxies is much more concentrated 
than that of the dissipationless dark matter particles. Sim¬ 
ulations predict that the mass loss of infalling subhaloes 


depends inversely on their halo-centric radius (e.g. 

Springel 

let al.|200111De Lucia et al.|200411Gao et al.|200411Xi 

e & Gao 
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2 Ran Li et. al 


2015). Thus, the halo mass to stellar mass ratio of satellite group catalog of Yang et al. (20071 which is constructed 


galaxies should increase as a function of halo-centric radius. 
Mapping the mass function of subhaloes from observations 
can provide important constraints on this galaxy evolution 
model. 

In observations, dark matter distributions are best mea¬ 
sured with gravitational lensing. For dark matter subhalos, 
however, such observations are challenging due to their rela¬ 
tive low mass compared to that of the host dark matter halo. 
The presence of subhalos can cause flux-ratio anomalies in 


multiply imaged lensing systems (Mao & Schneider 
Metcalf^^_^adn^[2001[ |Mao et al.||2004| |Xu et al. 


jan 


1998 


2009 


Nierenberg et al.|2014l, perturb the locations, and change 


the image numbers (Kneib et al. 1996 Kneib & Natara- 


2011), and disturb the surface brightness of extended 


Einstein ring/arcs (Koopmans 2005 Vegetti & Koopmans 
2009a|b Vegetti et al.|20i0 2012 l. However, due to the lim¬ 


ited number of high quality images and the rareness of strong 
lensing systems, only a few subhalos have been detected 
through strong lensing observation. Besides, strong lensing 
effects can only probe the central regions of dark matter 


haloes (Kneib & Natarajan 20111. Therefore, through strong 


lensing alone, it is difficult to draw a comprehensive picture 
of the co-evolution between subhalos and galaxies. 

Subhalos can also be detected in individual clusters 
through weak gravitational lensing or weak lensing combin- 


et al. 


ing strong lensing (e.g. Natarajan et al. 2007 2009 Limousin 


2005 2007 Okabe et al. 2014). In Natarajan et al. 


(20091, Hubble Space Telescope images were used to inves¬ 


tigate the subhalo masses of L* galaxies in the massive, lens¬ 
ing cluster, C10024+16 at z = 0.39, and to study the subhalo 


mass as a function of halo-centric radius. Okabe et al. (20141 


investigated subhalos in the very nearby Coma cluster with 
imaging from the Subaru telescope. The deep imaging and 
the large apparent size of the cluster allowed them to mea¬ 
sure the masses of subhalos selected by shear alone. They 
found 32 subhalos in the Coma cluster and measured their 
mass function. However, this kind of study requires very high 
quality images of massive nearby clusters, making it hard to 
extend such studies to large numbers of clusters. 

A promising alternative way to investigate the satellite- 
subhalo relation is through a stacking analysis of galaxy- 
galaxy lensing with large surveys. Different methods have 


been proposed in previous studies (e.g. Yang et al. 2006 


Li et al.||2013 Pastor Mira et ^|2011 Shirasaki||2015 ). Al 

though the tangential shear generated by a single subhalo is 
small, by stacking thousands of satellite galaxies the statisti¬ 
cal noise can be suppressed and the mean projected density 
profile around satellite galaxies can be measured. |Li et al(] 
(2014) selected satellite galaxies in the SDSS group catalog 


from (Yang et al. 20071 and measured the weak lensing signal 


around these satellites with a lensing source catalog derived 


using SDSS spectroscopic galaxies, the redMaPPer catalog 
relies on photometric cluster detections, allowing it to go to 
higher redshifts. As a result, there are more massive clusters 
in the redMaPPer cluster catalog. Therefore, we expect to 
signal-to-noise of the satellite galaxy lensing signals to be 
higher, enabling us to derive better constraints on subhalo 
properties. 

The paper is organized as follows. In section]^ we de¬ 
scribe the lens and source catalogs. In sectionwe present 
our lens model. In section]^ we show our observational re¬ 
sults and our best fit lens model. The discussions and con¬ 
clusions are presented in section Throughout the paper, 
we adopt a ACDM cosmology with parameters given by the 
WMAP-7-year data ( Komatsu et al. 12010 ) with Q,l = 0.728, 

= 0.272 and h = 7/o/(100kms“^Mpc“^) = 0.73. In this 
paper, stellar mass is estimated assuming a Chabrier (2003) 
IMF. 


2 OBSERVATIONAL DATA 

2.1 Lens selection and stellar masses 

We use satellite galaxies in the redMaPPer clusters as lenses. 
The redMaPPer cluster catalog is extracted from photomet¬ 


ric galaxy samples of the SDSS Data Release 8 (DR8, Aihara 


et al. 2011) using the red-sequence Matched-filter Proba¬ 
bilistic Percolation cluster Ending algorithm ( [Rykoff et al.| 
2014). The redMaPPer algorithm uses the 5-band (ugriz) 
magnitudes of galaxies with a magnitude cut i < 21.0 over 
a total area of 10,000 deg^ to photometrically detect galaxy 
clusters. 

redMaPPer uses a multi-color richness estimator A, de- 
Hned to be the sum of the membership probabilities over all 
galaxies. In this work, we use clusters with richness A > 20 
and photometric redshift Zduster < 0.5. In the overlapping 
region with the CFHT Stripe-82 Survey, we have a total of 
634 clusters. For each redMaPPer cluster, member galaxies 
are identified according to their photometric redshift, color 
and their cluster-centric distance. To reduce the contami¬ 
nation induced by fake member galaxies, we only use satel¬ 
lite galaxies with membership probability Pmem > 0.8. The 
redMaPPer cluster finder identifies 5 central galaxy candi¬ 
dates for each cluster, each with an estimate of the proba¬ 
bility Pcen that the galaxy in question is the central galaxy 
of the cluster. We remove all central galaxy candidates from 
our lens sample. For more details about redMaPPer cluster 


catalog, we refer the readers to Rykoff et al. (20141; Rozo & 


Rykoff (2014). 


Stellar masses are estimated for member galaxies in the 
redMaPPer catalog using the Bayesian spectral energy dis- 


from the CFHT Stripe82 Survey (Comparat et al. 20131. tribution (SED) modeling code iSEDfit (Moustakas et al. 


This was the first galaxy-galaxy lensing measurement of sub¬ 
halo masses in galaxy groups. However, the uncertainties of 
the measured subhalo masses were too large to investigate 
the satellite-subhalo relation as a function of halo-centric 
radius. 

In this paper, we apply the same method as |Li et al(] 
(2013 2014) to measure the galaxy-galaxy lensing signal for 
satellite galaxies in the SDSS redMaPPer cluster catalog 


(Rykoff et al. 2014 Rozo & Rykoff 2014). Unlike the the 


2013). Briefly, iSEDfit determines the posterior probability 
distribution of the stellar mass of each object by marginaliz¬ 
ing over the star formation history, stellar metallicity, dust 
content, and other physical parameters which influence the 
observed optical/near-infrared SED. The input data for each 
galaxy includes: the SDSS ugriz model fluxes scaled to the 
r-band cmodel flux; the 3.4- and 4.6-p.m “forced” WISE 
( Wright et al.j2010| photometry from Lang et al. (2014); 
and the spectroscopic or photometric redshift for each ob- 
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3 



Figure 1. The Mstar-zi distribution of lens galaxies, where 2 ; is 
the photometric redshift, and Mstar is in units of Mq. 


ject inferred from redMaPPer. We adopt delayed, exponen¬ 
tially declining star formation histories with random bursts 
of star formation snperposed, the Flexible Stellar Popula¬ 


tion Synthesis (FSPS) model predictions of Conroy et al. 
12009 [2010 1 , and the [Chabrier (20031 initial mass func¬ 
tion (IMF) from 0.1-100 Mq. For reference, adopting the 


Salpeter (19551 IMF would yield stellar masses which are 


on average 0.25 dex (a factor of 1.8) larger. We apply a stel¬ 
lar mass cut of Mstar > lO'^°h~'^M 0 to our satellite galaxy 
sample. In Fig[^ we show the Matar-^i distribution for our 
lens samples, where Zi is the photometric redshift of the 
satellite galaxy assigned by the redMaPPer algorithm. The 
low stellar mass satellite galaxies are incomplete at higher 
redshift, but they will not affect the conclusion of this paper. 


2.2 The Source Catalog 


The source catalog is measured from the Canada-France- 
Hawaii Telescope Stripe 82 Survey (CS82), which is an 
r—band imaging survey covering the SDSS Stripe82 region. 
With excellent seeing conditions — FWHM between 0.4 to 
0.8 arcsec — the CS82 survey reaches a depth of iab ~ 24.0. 
The survey contains a total of 173 tiles, 165 of which from 


CS82 observations and 8 from CFHT-LS Wide (Erben et al. 


2013). The CS82 fields were observed in four dithered ob¬ 


servations with 410s exposure. The 5a limiting magnitude 
is iAB ~ 24.0 in a 2 arcsec diameter aperture. 

Each CFHTLenS science image is supplemented by a 
mask, indicating regions within which accurate photome¬ 
try/shape measurements of faint sources cannot be per¬ 
formed, e.g. due to extended haloes from bright stars. Ac¬ 
cording to [Erben et al.| ( |2013| ), most of the science analysis 
are safe with sources with MASK^ 1. After applying all 
the necessary masks and removing overlapping regions, the 
effective survey area reduces from 173 deg^ to 129.2 deg^. 
We also require that source galaxies to have FITCLASS= 0, 
where FITCLASS is the flag that describes the star/galaxy 
classification. 

Source galaxy shapes are measured with the lensfit 
method ( Miller et al.|200T 20131, closely following the pro¬ 
cedure in Erben et al. (2009 20131. The shear calibration 


and systematics of the lensfit pipeline are described in detail 


in Heymans et al. (2012). The specific procedures that are 
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applied to the CS82 imaging are described in Erben et al. 
(2015, in preparation). 

Since the CS82 survey only provides the i-band images, 
the CS82 collaboration derived source photometric redshifts 


using the ugriz multi-color data from the SDSS co-add (An- 
nis et al.|2014 1, which reaches roughly 2 magnitudes deeper 


than the single epoch SDSS imaging. The photometric red- 
shift (photo-z) of the background galaxies were estimated 
using a Bayesian photo -2 code (BPZ, |Benftez||2000 Bundy 


et al.j 2015). The effective weighted source galaxy number 


density is 4.5 per arcmin . Detailed systematic tests for this 
weak lensing catalog are described in Leauthaud et al. 2015 
(in prep). 


2.3 Lensing Signal Computation 


In a galaxy-galaxy lensing analysis, the excess surface mass 
density , AE is inferred by measuring the tangential shear, 
7 t( 7 Z), where 


AE(i?) = 7 t(i?)E„it = E(< R) - E(i?), 


( 1 ) 


where E(< R) is the mean surface mass density within R, 
and E(i?) is the average surface density at the projected 
radius R, and Ecrit is the lensing critical surface density 


Ecrit — 


( 2 ) 


Ds 

4nG DiDis ’ 

where Dig is the angular diameter distance between the lens 
and the source, and Di and Dg are the angular diameter 
distances from the observer to the lens and to the source, 
respectively. 

We select satellite galaxies as lenses and stack lens- 
source pairs in physical radial distance R bins from 0.04 to 
1.5h“^Mpc. To avoid contamination from foreground galax¬ 
ies, we remove lens-source pairs with Zg — zi < 0.1, where Zg 
and zi are lens redshift and source redshift respectively. We 
have also tested the robustness of our results by varying the 
selection criteria for source galaxies. We find that selecting 
lens-source pairs with Zg — zi > 0.05 or Zg — zi > 0.15 only 
changes the final lensing signal by less than 7%, well below 
our final errors. 

For a given set of lenses, AE(J?) is estimated using 

^2l^2g ^crit(2:p 2s) 


AE(i?) = 




where 


Wls — ^crit 5 -^s) 


(3) 


(4) 


and Wn is a weight factor defined by Eq. ( 8 ) in Miller et al. 


(2013). w„ is introduced to account for the intrinsic distri¬ 


bution of ellipticity and shape measurement uncertainties. 

In the lensfit pipeline, a calibration factor for the mul¬ 
tiplicative error m is estimated for each galaxy based on 
its signal-to-noise ratio, and the size of the galaxy. Follow¬ 
ing Miller et al. (20131, we account for these multiplicative 


errors in the stacked lensing by the correction factor 


1 4- K{R) = 


EiEs'“^i«(4 + ms) 

Ei Es 


The corrected lensing measurement is as: 

led AE(i?) 


AE(i?)" 


1 -b K{R) 


(5) 


( 6 ) 
























































4 Ran Li et. al 


In this paper, we stack the tangential shears around 
satellite galaxies binned according to their projected halo¬ 
centric distance rp, and fit the galaxy-galaxy lensing signal 
to obtain the subhalo mass of the satellite galaxies. We de¬ 
scribe our theoretical lens models below. 


3 THE LENS MODEL 


The surface density around a lens galaxy E(i?) can be writ¬ 
ten as: 

E{R) = dx ; (7) 

and the mean surface density within radius R is 

E{<R) = ^ l\{u)udu, ( 8 ) 

where Pg,m is the density profile around the lens, and x 
is the physical distance along the line of sight. The excess 
surface density around a satellite galaxy is composed of three 
components: 


AS(i?) — ASaub(i?) + AEhost rp) -|- AEstar(I?) I (9) 


where AEaub(i?) is the contribution of the subhalo in which 
the satellite galaxy resides, AEhoat(I?, Cp) is the contribution 
from the host halo of the cluster/group, where rp is the 
projected distance from the satellite galaxy to the center of 
the host halo, and AEstar(7?) is the contribution from the 
stellar component of the satellite galaxy. We neglect the two- 
halo term, which is the contribution from other haloes on the 
line-of-sight, because this contr ibution is only im portant at 
R > 3h“^Mpc for clusters (see Shan et al.||2015 |. At small 
scales where the subhalo term dominates, the contribution 
of the 2-halo term is at least an order of magnitude smaller 
than the subhalo term (Li et al.|2009 1. 


3.1 Host halo contribution 


We assume that host haloes are centered on the central 
galaxies of clusters, with a density profile following the NFW 
( [Navarro et al. jpfft formula: 

p 0 ,host 


Phost (^) — ' 


( 10 ) 


(1 -I- r/ra,host)^(r/ra,host) ’ 

where rs^host is the characteristic scale of the halo and po.host 
is a normalization factor. Given the mass of a dark matter 
halo, its profile then only depends on the concentration pa¬ 
rameter c = r 2 oo/r’s,host, where r 2 oo is a radius within which 
the average density of the halo equals to 200 times the uni¬ 
verse critical density, pcrit. The halo mass M 200 is defined as 
M 200 = 47r/3r2oo(200pcrit). 

Various fitting formulae for mass-concentration rela¬ 
tions have been derived from N-body numerical simulations 


(e.g. Bullock et al.|2001 

Zhao et al. 

2003| Dolag et al.|2004| 

Maccio et al.|20071 [Zhao et af.|2009 

[Duffy et af.|2008||Neto 

et al. 

2007 

I. These studies find that the concentration de- 


creases with increasing halo mass. Weak lensing observations 
also measure this trend, but the measured amplitude of the 
mass-concentration relation is slightly smaller than that in 
the simulation (e.g. Mandelbaum et al.|2008 Miyatake et al. 


M^jShan et al.|2015| ). Since there is almost no degeneracy 

between the subhalo mass and the concentration (Li et al. 


2013 George et al. 20121, we expect that the exact choice 


of the mass-concentration relation should not have a large 
impact on our conclusions. Throughout this paper, we adopt 


the mass concentration relation from Neto et al. (20071: 


c = 4.67(M2oo/lO^'‘h"^M0 


( 11 ) 


We stack satellite galaxies with different halo-centric 
distance rp, and in host halos with different mass. Thus 
the lensing contribution from the host halo is an average 
of AEhost over rp and host halo mass M 200 . The host halo 
profile is modeled as follows. 

For each cluster, we can estimate its mass via the 


richness-mass relation of Rykoff et al. (20121: 


In 


M 200 


h^Q IO^^A/q 


= 1.48-f 1.06 ln(A/60) (12) 


In the redMaPPer catalog, each cluster has five possible cen¬ 
tral galaxies, each with probability Pcen. We assume that the 
average AE contribution from the host halo to a satellite can 
be written as: 


AEhost(P) = Ao AEhost(P|rpj, M2 qo)Pc( 


(13) 


where rpj is the projected distance between the satellite 
and the jth candidate of the host cluster center, Pcen,j is the 
probability of jth candidate to be the true center, and Nas.t 
is the number of stacked satellite galaxies. Aq is the only free 
parameter of the host halo contribution model. It describes 
an adjustment of the lensing amplitude. If the richness-mass 
relation is perfect, the best-fit Aq should be close to unity. 
Note that, the subhalo mass determination is robust against 
the variation of the normalization in richness-mass relation. 
If we decrease the normalization in Eq. [^by 20%, the best- 
fit subhalo mass changes only by 0.01 dex, which is at least 
15 times smaller than the Icr uncertainties of Msub (see table 
§. 


3.2 Satellite contribution 


In numerical simulations, subhalo density profiles are found 


to be truncated in the outskirts (e.g. Hayashi et al. 2003 


Springel et al.|[2008 Gao et aL||2004| jXie &: Gao||2015 1. In 

this work, we use a truncated NFW profile ( jBaltz et ah 
2009 Oguri fc Hamana|2011 tNFW, hereafter) to describe 
the subhalo mass distribution: 


Psub(r) — 


Po 


r/r^l -br/r,)2 -fr 


rt 


(14) 


where rt is the truncation radius of the subhalo, Vs is the 
characteristic radius of the tNFW profile and po is the nor¬ 
malization. The enclosed mass with x = r/vs can be written 
as: 


M{< x) = 47rpo’'s 


' 2(r2 -|- 1)2(1 -I- x){t^ + x^) 

X + l)x {x{x -|- 1) — R{x — 1)(2 -I- 3a;) — 2r‘*} 
-|-r(a; -I- l)('r^ -|- x^) {2(3r^ — l)arctan(a;/r) 

-l-r(r^ — 3) ln(r^(l + x)‘^l{R + 2 ;^))}], 


(15) 
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where t = rtjrs- We define the subhalo mass M^ub to be the 
total mass within rt- Given Msuh, Ts and rt, the tangential 
shear 7 t can be derived analytically (see the appendix in 
Oguri & Hamana (20111). 


Previous studies have sometimes used instead the 
pseudo-isothermal elliptical mass distribution (PIEMD) 
model derived by Kassiola & Kovner (1993l) for modeling 
the mass distribution around galaxies (e.g. Limousin et al. 


2009 Natarajan et al.|2009 Kneib fc Natarajan|2011 1. The 

surface density of the PIEMD model can be written as: 


E(i?) = 


l-Ro/Rt \-^Rl + R^ y/R'^t + R? 


(16) 


where Rq is the core radius, Rt is the truncation radius, 
and Eq is a characteristic surface density. The subhalo mass 
Msub, which is defined to be the enclosed mass with Rt, can 
be written as: 


yjRl-R^- ^R^t - + {Rt - Ro) 

(17) 

In this paper, we will fit the data with both tNFW and 
PIEMD models. 

Finally, the lensing contributed from stellar component 
is usually modeled as a point mass: 

AE.tar(i?) = ^ , (18) 

where Matar is the total stellar mass of the galaxy. 


_ 2'KT,oRoRt 

Rt — Ro 


4 RESULTS 

4.1 Dependence on the projected halo-centric 
radius 


To derive the subhalo parameters, we calculate the x de¬ 
fined as: 

x" = - AE(R0°'’“)(^')(AE(R,-) - AEiRjY’’^), 

^3 


(19) 

where AT,{Ri) and A'E{Ri)°’’‘‘ are the model and the ob¬ 
served lensing signal in the ith radial bin. The matrix Cij is 
the covariance matrix of data error between different radius 
bins. Even if the ellipticity of different sources are indepen¬ 
dent, the cross term of the covariance matrix could still be 
non-zero, due to some source galaxies are used more than 


once(e.g. Han et al. 2015 I. The covariance matrix can be rea¬ 
sonably calculated with bootstrap method using the survey 
data themselves! Mandelbaum et al.|2005 (. In the paper, We 
divide the CS82 survey area into 45 equal area sub-regions. 
We then generate 3000 bootstrap samples by resampling the 
45 sub-regions of CS82 observation data sets and calculate 
the covariance matrix using the bootstrap sample. Thus, the 
likelihood function can be given as: 


L oc exp • ( 20 ) 

We select satellite galaxies as described in section |2.1| 
and measure the stacked lensing signal around satellites in 
three Vp bins: 0.1 < < 0.3; 0.3 < rp < 0.6 and 0.6 < rp < 

0.9(in unit of h“^Mpc). 

For each bin, we fit the stacked lensing signal with 


Table 1. Flat priors for model parameters. M^ub is in units of 
and distances are in units of h“^Mpc. 



lower bound 

Upper bound 

Ao 

0.3 

2 

log Msub 

9 

13 

rt (tNFW) 

0 

0.35 

rs (tNFW) 

0 

0.06 

Ro (PIEMD) 

0 

0.02 

Rt (PIEMD) 

0 

0.35 


the 

the 


a Monte-Carlo-Markov-Chain(MCMC) method. For 
tNFW subhalo case, we have 4 free parameters: Msub 
subhalo mass; Vs, the tNFW profile scale radius; rt, the 
tNFW truncation radius, and Ao, the normalization factor 
of the host halo lensing contribution. 

We adopt flat priors with broad boundaries for these 
model parameters. We set the upper boundaries for rt and 
Ts to be the value of the viral radius and the scale radius 
of an NFW halo of 10^®h“^MQ. We choose these values be¬ 
cause the subhalo masses of satellite galaxy in clusters are 
expected to be much smaller than lO^^h^'^M© (Gao et al. 


20121 . 


For the PIEMD case, we also have 4 free parameters: 
Msub, Ro, Rt and Aq. We set the upper boundary of Rt 
to be the same as rt in the tNFW case. We set the upper 
boundary of Ro to be 20 kpc/h, which is much higher than 
the typical size of Ro in observations! Limousin et al.||2005 


Natarajan et al. 2009|). We believe our choice of priors is 


very conservative. The detailed choices of priors are listed 
in table [T] 

In Fig[^ we show the observed galaxy-galaxy lensing 
signal. Red dots with error bars represent the observational 
data. The vertical error bars show 1 cr errors estimated with 
the bootstrap resampling the lens galaxies. Horizontal error 
bars show the range of each radial bin. The lensing signals 
show clearly the characteristic shape described in figure 2 of 


Li et al. (20131. The lensing signal from the subhalo term 
dominates the central part. On the other hand, the con¬ 
tribution from the host halo is nearly zero on small scale, 
and decreases to negative values at intermediate scale. This 
is because Ehoat(R) becomes increasingly larger compared 
to Ehost(< R) at intermediate R. At radii where R > rp, 
AEhoat(R) increases with R again. At large scales where 
R » rp, AT,{R,rp) approaches AEhoat(R, 0). 

The solid lines show the best-fit results with the tNFW 
model. Dashed lines of different colors show the contribution 
from different components. The best-fit model parameters 
are listed in table Note that, the value of the first point 
in the [0.6,0.9] rp bin is very low, which may be due to the 
relatively small number of source galaxies in the inner most 
radial bin. We exclude this point when deriving our best fit 
model. For comparison, we also show the best-fit parameters 
including this point in table 

Fig.|^shows an example of the MCMC posterior distri- 
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bution of the tNFW model parameters for satellites in the 
[0.3, 0.6]h~^Mpc rp bin. The constraints on the subhalo mass 
logMsub are tighter than that in Li et al. (2014 1 (~ ±0.7). 
The amplitude of Aq is slightly smaller than unity, imply¬ 
ing that the clusters may be less massive than predicted 
by the mass-richness relation. However, we are still unable 
to obtain significant constraints on the structure parame¬ 
ters of the sub-halos. In principle, the density profile cut¬ 
off caused by tidal effects can be measured with tangen¬ 
tial shear. However, the galaxy-galaxy lensing measurement 
stacks many satellites, leading to smearing of the signal. 
With the data used here, the tidal radius as a free parameter 
are not constrained. Some galaxy-galaxy lensing investiga¬ 
tions introduced additional constraints in order to estimate 
the tidal radius. For example, Gillis et al. (20131 assumed 
that galaxies of the same stellar mass but in different envi¬ 
ronments have similar sub-halo density profiles except the 
cut-off radius. With this additional assumption, they ob¬ 
tained rtidai /»'200 = 0.26 ± 0.14. During the review process 
of our paper, Sifon et al. (20151 posted a similar galaxy- 
galaxy lensing measurement of satellite galaxies using the 
KiDs survey, and they also found that their data cannot 
distinguish models with or without tidal truncation. 

In Fig[4l we show the fitting results of the PIEMD 
model. The best-fit parameters are listed in table For 
reference, we also over-plot the theoretical prediction of the 
best-fit tNFW model with blue dashed lines. Both models 
provide a good description of the data. The best-fitted Maub 
and Aq of the two models agree well with each other, show¬ 
ing the validity of our results. 

In numerical simulations, subhalos that are close to host 


halo centers are subject to strong mass stripping(Springel 
et al.||2001[ |De Lucia et al.||2004| |Gao et al.||2004| |Xie fc 


Gao|2015 1. The mass loss fraction of subhalos increases from 


30% at r 2 oo to 90% at 0 . 1 r 2 oo (IGao et al. 1 120041 IXie 


fc Gao|[2015[). From galaxy-galaxy lensing in this work, we 


also find that the Maub of [0.6, 0.9] Vp bin is much larger 
than that in [0.1,0.3] Vp bin (by a factor of 18). In Figj^ 
we plot the subhalo mass-to-stellar mass ratio for three rp 
bins. The Maub/Mstar ratio in the [0.6,0.9]h“^Mpc rp bin is 
about 12 times larger than that of the [0.1,0.3]h“^Mpc bin. 
If we include the first point in the [0.6, 0.9]h“^Mpc bin, the 
Maub/Mstar of the tNFW model decreases by 40%. 

For reference, we over-plot the Maub/Matar-^p predicted 
by the semi-analytical model of Guo et al. ( 20111). We adopt 


the mock galaxy catalog generated with the [Guo et al. 


(20111 model using the Millennium simulation (Springel 


et al. 


_ 2006 1. We select mock satellite galaxies with stel¬ 
lar masses Matar > lO^'^h^^M© from clusters with M200 > 
lO^^/h~^M 0 . The median Maub/Matar-rp relation is shown 
in Figj^ with a black solid line. The green shaded region 
represents the parameter space where 68 % of mock satellites 
distributes. We see the semi-analytical model predicts an in¬ 
creasing Maub/Matar with Tp, but with a flatter slope than in 
our observations. Note that we have not attempted to recre¬ 
ate our detailed observational procedure here, so source and 
cluster selection might potentially explain this discrepancy. 
Particularly relevant here is the fact that our analysis relies 
on a probability cut Pmem >0.8 for satellite galaxies, which 
implies that ~ 10 % of our satellite galaxies may not be true 
members of the clusters, but galaxies on the line of sight. 
This contamination is difficult to eliminate completely in 



i?/(h^^Mpc) 


Figure 2. Observed galaxy-galaxy lensing signal for satellite 
galaxies as a function of the radius. The top, middle and bot¬ 
tom panels show results for satellites with rp in the range 
[0.1, 0.3]h“lMpc, [0.3, 0.6]h“^Mpc and [0.6, 0.9]h“^Mpc, respec¬ 
tively. Red points with error bars represent the observational 
data. The vertical error bars show the 1 a bootstrap error. The 
horizontal error bars show the bin size. Black solid lines show the 
best-fit tNFW model prediction. Dashed lines of different colors 
show the contribution from the subhalo, the host halo, and the 
stellar mass respectively. 


galaxy-galaxy lensing, because the uncertainties in the line- 
of-sight distances are usually larger than the sizes of the 
clusters themselves. In Li et al. (20131, we used mock cata¬ 
logs constructed from N^-body simulations to investigate the 
impact of interlopers. It is found that 10% of the galaxies 
identihed as satellites are interlopers, and this introduces 
a contamination of 15% in the lensing signal. We expect 
a comparable level of bias here, as shown in Appendix A. 
It should be noted, however, that the average membership 
probability of our satellite sample does not change signif¬ 
icantly with rp, implying that the contamination by fake 
member galaxies is similar at different rp. We therefore ex¬ 
pect that the contamination by interlopers will not lead to 
any qualitative changes in the shapes of the density profiles. 
A detailed comparison between the observation and simula¬ 
tion, taking into account the impact of interlopers, will be 
carried out in a forthcoming paper. 
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Figure 3. The contours show 68% and 95% confidence intervals for the tNFW model parameters, Mgub, r^, rj and Aq. Results are 
shown for satellites with rp = [0.3, 0.6]h“^Mpc. The last panel of each row shows the ID marginalized posterior distributions. Note that 
the plotting range for and rt is exactly the same as the limits of our prior, so we actually do not have much constraints for these two 
parameters, except that high values are slightly disfavored for both. 


4.2 Dependence on satellite stellar mass 


In the CDM structure formation scenario, satellite galaxies 
with larger stellar mass tend to occupy larger haloes at infall 



satellite galaxies with larger stellar mass reside in more mas¬ 
sive subhalos. 

To test this prediction, we select all galaxies with rp in 
the range [0.3, 0.9] h“^Mpc, and split the them into two 
subsamples: 10 < log(Mstar/h“^MQ) < 10.5 and 11 < 
log(Mstar/h“'^M 0 ) < 12. The galaxy-galaxy lensing signals 
of the two satellite samples are shown in Fig|^ It is clear that 
at small scales, where subhalos dominate, the lensing signals 
are larger around the more massive satellites. The best-fit 
subhalo mass for the high mass and low mass satellites are 
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Table 2. The best-fit values of the tNFW model parameters for the stacked satellite lensing signal in different Vp bins. logMsub? Aq 
are the best-fit values for the subhalo mass, and the normalization factor. We do not show the best-fit of Vs and rt as they are poorly 
constrained. Nsat is the number of satellite galaxies in each bin. (logMstar) is the average stellar mass of satellites. All errors indicate 
the 68% confidence intervals. Masses are in units of h~^MQ. In our fiducial fitting process, we exclude the first point in [0.6, 0.9]h“^Mpc 
Tp bin as a outlier (see Fig[§- For comparison, the bottom row of the table shows the fitting results including this first point. 


Rp range 

(log(Matar)> 

log Mgub 

^0 

-^sub/-^star 

A^sat 

i^i) 

[0.1, 0.3) 

10.68 

1 -| 07+0-35 

O 

bo 

o 

1 + 

op 

b b 

4 4^+0 03 

3963 

0.33 

[0.3, 0.6) 

10.72 

ii.92i®:« 

0.86l°;“ 

1 7 00+6.98 

1 i •^o_g 34 

2507 

0.29 

[0.6, 0.9) 

10.78 

12.64t“;[? 

O O 

bb 
+ 1 

00 

o 

<0.4U_;L9.09 

577 

0.24 

[0.6, 0.9)* 

10.78 

12.491®;)® 

o o 
bb 
+ 1 

00 

o 

54 «4+16.58 

04.D^_;L5.80 

577 

0.24 


Table 3. best-fitted parameters with PIEMD model. 


Vp range 

(log(Mstar)} 

log Msuh 

Ao 

A^sub / Afstar 

[0.1, 0.3) 

10.68 

1 1 oq+0.55 

li.CiU_Q 57 

0.80l°o®) 

0-'^-5.59 

[0.3, 0.6) 

10.72 

11 76+0.33 
li. /D_q 32 

0.861?-®) 

14.401?;?? 

[0.6, 0.9) 

10.78 

12.80lg;)? 

0-791?:?) 

110.981??;!? 

[0.6,0.91* 

10.78 

12.84l?;« 

0-791?:?) 

119.741??;?? 


log(M,,b/h-iM0) = 11.141“;?® = 19.5ll?'j) 

and log(Msub/h“^M 0 ) = 12.381qi? (Msub/Matar = 

2I.II77 respectively. 


5 DISCUSSION AND CONCLUSION 

In this paper, we present measurements of the galaxy-galaxy 
lensing signal of satellite galaxies in redMaPPer clusters. We 
select satellite galaxies from massive clusters (richness A > 
20) in the redMaPPer catalog. We fit the galaxy-galaxy lens¬ 
ing signal around the satellites using a tNFW profile and a 
PIEMD prohle, and obtain the subhalo masses. We bin satel¬ 
lite galaxies according to their projected halo-centric dis¬ 
tance rp and find that the best-fit subhalo mass of satellite 


galaxies increases with rp. The best-£t log(Maub/h“^M 0 ) 
for satellites in the rp = [0.6, 0.9]h““Mpc bin is larger than 
in the Vp = [0.1,0.3]h~'^Mpc bin by a factor of 18. The 
Afsub/Mstar ratio is also found to increase as a function of 
rp by a factor of 12. We also find that satellites with more 
stellar mass tend to populate more massive subhalos. Our 
results provide evidence for the tidal striping effects on the 
subhalos of red sequence satellite galaxies, as expected on 
the basis of the CDM hierarchical structure formation sce¬ 
nario. 

Many previous studies have tried to test this theoretical 
prediction using gravitational lensing. Most of these previ¬ 
ous studies focus on measuring subhalo mass in individual 
rich clusters. For example, Natarajan et al. (20091 report the 
measurement of dark matter subhalos as a function of pro¬ 
jected halo-centric radius for the cluster C10024-I-16. They 
found that the mass of dark matter subhalos hosting early 
type L* galaxies increases by a factor of 6 from a halo-centric 
radius r < 0.6 h““Mpc to r ~ 4 h“^Mpc. In recent work. 


Okabe et al. (20141 present the weak lensing measurement 


of 32 subhalos in the very nearby Coma cluster. They also 
found that the mass-to-light ratio of subhalos increase as a 
function of halo-centric radius: the M/Li/ of subhalos in¬ 
creases from 60 at 10' to about 500 at 70'. 

Our work is complementary to these results. In the pa¬ 
per, we measure the galaxy-galaxy lensing signal of subhalos 
in a statistical sample of rich clusters. Our results lead to 
similar conclusions and show evidence for the effects of tidal 
striping. However, due to the statistical noise of the current 
survey, there is still a large uncertainties in the measurement 
of Msub/Matar. 

Theoretically, the galaxy - halo connection can be es¬ 
tablished by studying how galaxies of different properties 
reside in dark matter halos of different masses, through 
models such as the halo occupation models Jing et al. (e.g. 


1998P; Peacock & Smith (e.g. 20001; Berlind & Weinberg 


(e.g. [ 2002 1 , the conditional luminosity functions models (e.g. 


© 0000 RAS, MNRAS 000, 000-000 






























9 


U 

a 

© 




W 

<1 



"02 ^4 06 08 ito 1^2 L4 



Figure 4. This figure is similar to Fi^ Solid lines represent the 
best-fit PIEMD model prediction and blue dashed lines are the 
predictions of the tNFW model. 



Figure 5. The subhalo mass-to-stellar mass ratio for galaxies in 
different Vp bins. Vertical error bars show the 68% confidence in¬ 
terval. Empty circles show the best-fitted Mgub/^star values for 
[0.6, 0.9]h“^Mpc Tp bin including the inner most observational 
point (see the bottom panel in Fig[§ . Horizontal error bars show 
the bin range of rp. The green shaded region represents the pa¬ 
rameter space where 68% of the semi-analytical satellite galaxies 
distributes. The semi-analytical model also predicts an increasing 
A^sub/^star with Tp, but with a flatter slope. 



Figure 6. Observed galaxy-galaxy lensing signal for satellite 
galaxies in different stellar mass bins at fixed . The legend shows 
the log(Mstar/h~^MQ) range for different data points. The solid 
lines are the best-fit tNFW models. 


Yang et al. 

2003 

1, (sub)halo-abundance-matching ( e.g. Vale 

& Ostriker 

2004 

IConroy et al.|2006||Wang et al.|20061|Yang 

et al. 2012 

Chaves-Montero et al. 2015 

1, and empirical mod- 


els of star formation and mass assembly of galaxies in dark 
matter halos (e.g. Lu et al. 2014 20151. In most of these 
models, the connection is made between galaxy luminosities 
(stellar masses) and halo masses before a galaxy becomes a 
satellite, i.e. before a galaxy and its halo have experienced 
halo-specific environmental effects. The tidal stripping ef¬ 
fects after a halo becomes a sub-halo can be followed in dark 
matter simulations. This, together with the halo-galaxy con¬ 
nection established through the various models, can be used 
to predict the subhalo - galaxy relation at the present day. 


During the revision of this paper, Han et al. (20161 posted a 


theoretical work on subhalo spatial and mass distributions. 
In Sec. 6.4 of their paper, they presented how our lensing 
measurements can be derived theoretically from a subhalo 
abundance matching point of view. With future data, the 
galaxy-galaxy lensing measurements of subhalos associated 
with satellites are expected to provide important constraints 
on galaxy formation and evolution in dark matter halos. 

The results of this paper also demonstrate the promise 


of next generation weak lensing surveys. In Li et al. (20131, 


it is shown that the constraints on subhalo structures and 
ATjub/Mstar can be improved dramatically with next gener¬ 
ation galaxy surveys such as LSST, due to the increase in 
both sky area (17000 deg^) and the depth (40 gal/arcmin^), 
which is crucial in constraining the co-evolution of satellite 
galaxies and subhalos. The space mission, such as Euclid, 
will survey a similar area of the sky (20000 deg^) but with 
much better image resolutions (FWHM 0.1” versus 0.7” for 
LSST), which is expected to provide even better measure- 
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ments of galaxy-galaxy leasing. The method described here 
can readily be extended to these futnre surveys. 
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APPENDIX A: FAKE MEMBER 
CONTAMINATION 


The uncertainties of photometric redshifts cause the error 
in lens selection. The line-of-sight field galaxies may be 
misidentified as cluster members. These held galaxies are 
free of tidal disruption, causing an overestimation of the 
Msub/Mstar for Satellite galaxies. Rozo et al. (20151 investi¬ 
gate the cluster member identihcation, and hnd that the 
photometric membership estimator agrees with the spec¬ 
troscopic membership with 1% precision. In this work, 
we use satellite galaxy with membership probability larger 
than 80%, which give us a mean probability of satellite 
~ 80% — 90%. Thus, we conclude that our satellite selec¬ 
tion may face a 10%-15% contamination from the light sight 
held galaxies. According to E et al.| ( |2013[ ), 10% interlop¬ 
ers contribute to about 15% of the total signal in the in¬ 
ner part. The accurate estimation of the contamination re¬ 
quires a detailed comparison with mock catalog constructed 
in proper simulations, which is not available at the moment. 
We thus made the following test to explore the effects of 
the interlopers. Instead of choosing satellite galaxies with 
Pmem > 80%, we choose satellite galaxies with Pmem > 60%, 
and re-calculate the lensing signal. We show the comparison 
of the lensing signal in Fig |Al[ The later Pmem selection will 
increase the interloper fraction from 10%-15% to 20%-25%. 
We hnd that the change in the lensing is only about 10% at 
the subhalo dominating region, smaller than the statistical 
error bars. 
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Figure Al. Observed galaxy-galaxy lensing signal for satellite 
galaxies with different Pmem selection. 


APPENDIX B: PHOTOMETRIC REDSHIFT 
BIAS TEST 


The biased photometric redshifts can cause a systematic er¬ 
ror in the measured lensing signal. We make the follow¬ 
ing test to explore this uncertainty. We split the source 
galaxy sample in two, the high-z {z > 0.7) and the low- 
z {z < 0.7) samples. We then compute AT,{R) for the 
same set of lenses with each half of the sources. We ht 
A'Ehigh-z = (1 -f 6)AEiow-z where b is a free parameter, 
which is used to estimate the bias introduced by photo-z 
error. If the source redshifts are free from bias, we will have 
b=0. In this test, the lens sample can be any galaxies as 
long as their redshift distribution is the same with that of 
our satellite samples in the paper. To maximum the lensing 
signal, we thus select galaxies from the LOWZ and GMASS 
galaxies from SDSS-BOSS surveyl Dawson et al.|20l3 ), and 
require them to have the same redshift distribution as the 
satellite galaxies in our paper. The results are shown in 
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Figure Bl. Observed galaxy-galaxy lensing signal with high-z 
sources (blue) and low-z sources (red). The lenses are selected 
from the LOWZ and CMASS galaxies. The lower panel shows 
the difference between results of the two sample of sources. 


Fig |Bl[ We find that the b parameter is consistent with 0 
within the statistical errors (6 = —0.05 ± 0.19). 
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